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Abstract 

The quark model does not exclude states composed of more than three quarks, 
like pentaquark systems. Controversial evidence for such states has been pub- 
lished in the last years, in particular: for a strange pentaquark 0(1540)"''; for 
a double-strange state, the S(1862) , subsequently called $(1860) ; and for 
a charmed state, the 0c (3 100)°. If confirmed, a full pentaquark family might 
exist; such pentaquark states could be produced in e^e~ annihilations near the 
Z energy. In this paper a search for pentaquarks is described using the DELPHI 
detector at LEP, characterized by powerful particle identification sub-systems 
crucial in the separation of the signal from the background for these states. At 
95% CL, upper limits are set on the production rates {N) of such particles and 
their charge-conjugate state per Z decay: 
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1 Introduction 

Pentaquark is a name given to describe a bound state of four quarks and one antiquark, 
e.g. uudds. The quark model docs not exclude such states. Several models predict the 
multiplet structure and characteristics of pentaquarks, for example the chiral soliton 
model, the uncorrelated and correlated quark models, the thermal model, lattice QCD 
etc. [1]. The current theoretical description of possible pentaquarks is very rich, but it 
does not provide a unique picture of the pentaquark characteristics. Furthermore, lattice 
calculations give very different predictions as to whether pentaquarks exist and, if they 
do, what mass and parity they have. 

Pentaquark states were first searched for in the 60's but the few, low statistics, pub- 
lished candidates were never confirmed [2]. More recent experimental evidence [3], how- 
ever, may suggest the existence of pentaquark systems. The first possible candidate is^ 
the 6(1540)+, with mass of (1.54 ± 0.01) GeV/c^ width smaller than 1 MeV/c^ and 
strangeness S=-|-l, consistent with being made of the quarks uudds. This evidence is still 
controversial as is that for the other pentaquark states discussed in this Letter (see [4] 
and references therein). 

Subsequently, evidence for another exotic baryon, doubly charged and with double 
strangeness, the S(1862) (subsequently called $(1860) , see [5]), has been claimed 
by the CERN experiment NA49 [6], with mass of (1862 ± 2) MeV/c^ 

Later, the DESY experiment HI has reported a signal for a charmed exotic baryon 
in the pD*~ channel [7], the 0c(31OO)'^. This resonance was reported to have a mass of 
(3099 ±3 (stat) ±5 (syst)) MeV/c^ and a measured width compatible with the experimen- 
tal resolution. It was interpreted as an anti-charmed baryon with a minimal constituent 
quark composition of uuddc. Several experiments tried to verify this finding [4]. The 
ZEUS collaboration for instance challenged the results of HI; even with a larger sample 
of D*+ mesons, such a narrow resonance was not observed [8] . 

Isospins and 1 are both possible for pentaquarks; isospin 1 would lead to three charge 
states ©°, ©+ and Thus the search is for a family of pentaquarks. 

This paper reports on the results of a search for pentaquark states in hadronic Z decays 
recorded by DELPHI. In a similar analysis, ALEPH [9] did not observe significant signals. 
The powerful particle identification characterizing the DELPHI detector might facilitate 
this search, since this feature helps in detecting and separating from the background some 
decay states of pentaquarks. 

The article is organised as follows. After a short description of the subdetectors used 
for the analysis (Section 2), Section 3 presents the results of a search for pentaquarks in 
the pK*^ (the ©"'") and the pK+ (the 0"'"+) channels. Section 4 presents a search for a 
doubly- charged, doubly-strange pentaquark (the $(1860) ). Section 5 presents a search 
for a charmed pentaquark (the 0c(31OO)°). A summary is given in Section 6. 

2 The Detector 

The DELPHI detector is described in detail in [10], and its performance is analysed 
in [11]. 

The present analysis relies mostly on information provided by the central tracking 
detectors and the Barrel Ring Imaging Cherenkov Counter (BRICH): 

^Charge conjugated states are implied throughout this paper. 
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• The micro Vertex Detector (VD) consists of three layers of sihcon strip detectors at 
radii^ of 6.3 cm, 9.0 cm and 10.9 cm. R(f) is measured in all three layers. The first 
and third layers also provide z information (from 1994 on). The 9 coverage for a 
particle passing all three layers is from 44° to 136°. The single point precision has 
been estimated from real data to be about 8 fim in R(j) and (for charged particles 
crossing perpendicular to the module) about 9 /xm in z. 

• The Inner Detector (ID) consists of an inner drift chamber with jet chamber geom- 
etry and 5 cylindrical MWPC (straw tube from 1995 on) layers. The jet chamber, 
between 12 and 23 cm in R and from 23° to 157° in 6 (15°-165° from 1995 on), 
consists of 24 azimuthal sectors, each providing up to 24 i?0 points. 

• The Time Projection Chamber (TPC) is the main tracking device. It provides up 
to 16 space points per particle trajectory for radii between 40 cm and 110 cm. The 
precision on the track elements is about 150 fim in R(j) and about 600 /xm in z. A 
measurement of the specific energy loss dE/dx of a track is provided with a resolution 
of about 6.5%, providing charged particle identification up to a momentum of about 

1 GeV/c. 

• The Outer Detector (OD) is a 4.7 m long set of 5 layers of drift tubes situated at 

2 m radius to the beam which provides precise spatial information in R(f). 

• The Barrel Ring Imaging Cherenkov Counter (BRICH) is the main DELPHI detector 
devoted to charged particle identification. It is subdivided into two halves (z > and 
2; < 0) and provides particle identification using Cherenkov radiation produced in a 
liquid or a gas radiator. This radiation, after appropriate focusing, is transformed 
into photoelectrons in a TPC-like drift structure and the Cherenkov angles of the 
track in both media are determined. The BRICH provides particle identification in 
the momentum range from 0.7 GeV/c to 45 GeV/c. 

The DELPHI tracking system was completed by two tracking chambers (FCA and FCB) 
in each forward region. 

To compute the selection efficiency of the various channels studied, Z — > gg events were 
simulated using the JETSET parton shower generator [12] and then processed through 
the DELPHI simulation program, DELSIM, which models the detector response. The 
simulated events passed through DELSIM were then processed by the same reconstruction 
program as used for the data, DELANA [11]. The amount of simulated events is more 
than twice the real data. 

For the 0+, 0++ and $(1860) searches, the data recorded during the LEPl oper- 
ation in the years 1991 to 1995 were used. For the 9c(3100)° search, the analysis was 
restricted to the years 1994 and 1995, the two highest luminosity years of LEPl, with all 
DELPHI particle identifiers fully operational. 

3 Search for Strange Pentaquarks in the pK system 

The state 0+ can be detected through its decay into pK° pairs; the state 0"*"+ could 
be detected through its decay into pK+. Therefore the invariant mass distributions of 
pK° and pK"^ pairs in hadronic Z decays were studied. These were compared with the 

pK^ s])ectrum. where the A(1520) is observed. 

■^In the standard DELPHI coordinate system, the z axis is along the electron beam direction, the x axis points towards 
the center of LEP, and the y axis points upwards. T he polar angle to the z axis is called 6 and the azimuthal angle around 
the z axis is called (/>; the radial coordinate is R= \/x'^ +y'^. 
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3.1 Event selection 

Hadronic Z decays for this analysis were selected by requiring at least four recon- 
structed charged particles and a total energy of these particles (assuming the pion mass) 
larger than 12% of the centre-of-mass (cm.) energy. The charged-particle tracks had 
to be longer than 30 cm, with a momentum larger than 400 MeV/c and a polar angle 
between 20° and 160°. The polar angle of the thrust axis, Othmst, was computed for each 
event and events were rejected if | cos9thrust\ was greater than 0.95. A total of 3.4 million 
hadronic events were selected. 

The selection efficiency for hadronic events was estimated using the simulation, and 
found to be larger than 95% within the angular acceptance. 

In order to search for the pentaquark states, the pK*^, pK^ and pK+ invariant mass 
spectra were constructed using identified particles. Particle identification was performed 
combining dE/dx and BRICH information. According to the quality of particle identifi- 
cation the tagging categories loose, standard and tight are distinguished for each particle 
species as well as for so-called "heavy" tag, which severely reduces the fraction of charged 
pions. To further improve the quality of particle identification for a track of given mo- 
mentum and (assumed) particle type it was required that information from the detectors 
specified in Table 1 was present. Only in the years 1994 and 1995 all particle identifi- 
cation detectors were fully operational; the identification was essentially only based on 
TPC during the years 1991 to 1993. 





momentum range in GeV/c 


0.3 - 0.7 


0.7-0.9 


0.9 - 1.3 


1.3- 2.7 


2.7-9.0 9.0-16.0 


16.0 - 45.0 


TT 


TPC 


LRICH S 


GRICH S 


K 


TPC 


LRICH S 


GRICH V 

+ 

LRICH S 


GRICH S 


P 


TPC 


TPC 

+ 

LRICH V 


LRICH S 


GRICH V 

+ 

LRICH S 


GRICH V 


GRICH S 



Table 1: Momentum ranges for particle identification: TPC denotes identification using 
the dE/dx measurement of the TPC, LRICH S (V) denotes identification using a signal 
(veto) of the hquid RICH, and correspondingly GRICH for the gas RICH. 

A particle was taken to be a proton if it was tightly tagged or fulfilled the standard 
tag by identification from ionization loss in the TPC. Kaons were required to be tightly 
tagged in the momentum ranges p < 3.5 GeV/c and p > 9.5 GeV/c. In the intermediate 
momentum range kaons were also identified by a tight heavy particle tag [13] combined 
with at least a standard kaon tag. 

3.1.1 Description of the invariant mass spectra 

In the present analysis, the mass spectra were described by a distribution function, 
/(M, a), of the invariant mass M. The parameters a were determined by a least squares 
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fit of tlie function to tlie data. Tlie function /(M, a) was composed of two parts: 

f{M,a)^f{M,a)+f{M,a), (1) 

corresponding to the signal and to the background respectively. The signal function, 
f^{M,a), described the resonance signals in the corresponding invariant mass distribu- 
tions. It has the form: 

f{M,a)^aixR{M,a2,a3), (2) 

where R is either a non-relativistic Breit-Wigner or a normalised Gaussian function ac- 
counting for the resonance production; a2 and 03 are respectively the fitted peak RMS 
width and mass m. The background term, f^{M,a), was taken to be a third order 
polynomial in M. 

3.2 Analysis of the pK^ channel 

The invariant mass distribution for pK*^ pairs was first studied. K° candidates were 
obtained from the fit of charged particle tracks of opposite charge consistent with the 
pion hypothesis, as described in [11]. The 7r"'"7r~ invariant mass is shown in Figure 1 (a). 

The pK° mass distribution is displayed in Figure 1 (b), for an invariant mass 
between 0.45 GeV/c^ and 0.55 GeV/c^. No signal is visible in the G"*" mass region; the 
simulation accounts very well for the data over the whole mass spectrum. 




Figure 1: (a) 7i~^7i invariant mass (b) pK° mass spectrum. The histogram represents 
the simulation, while the points represent the data. 

To set the hmit on the 0+ production, the fitting procedure as described above 
was applied, modeling a possible signal by a Gaussian function with a central value 
of 1.54 GeV/c^ and a RMS width''^ of 10 MeV/c^ equal to the resolution. 

The pK^ selection efficiency was estimated from a Monte Carlo generated sample of 
©"•" events to be (6.4 ± 0.3)%. The error is dominated by the systematic uncertainties 
coming from reconstruction and proton identification. 

'^Throughout the paper, if the fit is done by a Brcit-Wigiicr function the width indicates the value of the F parameter, 
while in the case of a Gaussian function it indicates the RMS error a. 
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The estimated number of events in the signal region was —20 ± 64 (stat). The corre- 
sponding upper Hmit, at 95% CL, on the average production rate per hadronic event of 
the e+ is: 

where the systematic uncertainty was added in quadrature to the statistical error. The 
result has been corrected for the branching fraction Br {Kg t:^-k~). 

3.3 Analysis of the pK~ and pK+ channels 

The search for a possible was made in the pK+ channel, after investigation of the 
channel pK~, where the presence of the A(1520) resonance allows the pK~(K"'") selection 
efficiency to be measured in the region of interest. Figure 2 (a) shows the pK~ invariant 
mass spectrum. A clear A(1520) signal is observed at the expected mass. It has been 
checked that there arc no prominent reflections from known particle decays in the pK~ 
mass spectrum. In addition pK^ combinations in which the combined with any 
identified K+ had a mass in the (j) (1020) region were discarded. The total excess in the 
A(1520) region, measured from the fit to the mass spectrum of Figure 2 (a) is of: 

(?^A(i520)) = 2 130 ± 450 events, (3) 

with a mass of 1.520 ± 0.002 GeV/c^ and a width of 0.010 ± 0.004 GeV/c^ compatible 
with the experimental resolution. The per degree of freedom is 1.4. The A(1520) 
selection efficiency determined from the simulation is (12.8 dz 0.5)%. This corresponds to 
an average A(1520) production rate per hadronic event of 0.0217 ± 0.0046 (stat) to be 
compared with the published value [5] of 0.0224 ± 0.0027. 

The invariant mass spectrum for pK+ pairs, obtained using the same cuts, is plotted 
in Figure 2 (b) . No significant peak is visible; the per degree of freedom of the fit to 
the background function only is 2.1. 




M(pK'), GeV/c' M(pK*), GeV/c' 

Figure 2: (a) Differential pK~ and (b) pK"*" mass spectra. The lines represent the fits 
described in the text. 

An upper limit for the average production rate of the 0"'""'" can be determined over 
the range of its mass estimates (1.45 to 1.65 GeV/c^), assuming the same efficiency as 
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for the A(1520). It should be taken into account that, since the A(1520) can decay into 
a charged pair and into a neutral pair as well, essentially with the same probability, the 
sensitivity to decay channels of the is twice that of the A(1520). 

A fit to the form (1) was performed by varying the mass between 1.45 GeV/c^ and 
1.65 GcV/c^ in steps of 5 MeV/c^, and by imposing a RMS width of 10 McV/c^ (the 
expected experimental resolution). Limits at 95% CL were then calculated as a function 
of the mass, yielding a maximum signal of 350 ± 187 events. The systematic uncertainties 
on the production rate of a G"*""*" with a mass close to the A(1520) mass, can be expected 
to be of the same order as those of the A(1520) production rate, which were estimated to 
be of 16% [14]. Such systematics are therefore negligible with respect to the error from 
the fit. A general limit 

{Ne++) < 1.6 X 10-^ 

for the mass region between 1.45 GeV/c^ and 1.65 GeV/c^ is obtained. This hmit is 
higher than what could be expected given the sensitivities, due to the about 2a statistical 
fiuctuations in the mass region between 1.52 GeV/c^ and 1.58 GeV/c^. 

4 Search for Doubly Charged and Doubly Strange 
Pentaquarks in the E~7v~ system 

The exotic baryons with double charge and double strangeness were searched through 
the decay into into S~7r~. The hadronic Z decays sample for this analysis is the same as 
described in section 3.1; it corresponds to a total of 3.4 million hadronic events after the 
cuts, recorded in the years 1991 to 1995. 

4.1 Reconstruction 

The S~ hyperon was reconstructed through the decay S~ — > Avr^. For this, all 
candidates, i.e., all pairs of oppositely charged particles, were considered as A candidates. 
For each pair, the higher momentum particle was assumed to be a proton and the other 
a pion, and a vertex fit performed using the standard DELPHI search algorithm [11]. 

The p7r~ invariant mass is shown in Figure 3 (a). The A candidates were selected 
by requiring an invariant mass M(p7r~) between 1.100 GeV/c^ and 1.135 GeV/c^, a 
probability of the vertex fit larger than 10^^ and a decay length from the interaction 
point greater than 0.2 cm in the plane transverse to the beam. 

A constrained multivertex fit was performed on each S~ candidate decaying into Att" 
[15]. The 16 measured variables in the fit were the five parameters of the helix parame- 
terization of each of the three charged particle tracks and the z coordinate of the beam 
interaction point (the x and y coordinates were so precisely measured that they could be 
taken as fixed). The fitted variables were the decay coordinates of the S~ and A. 

The fit constrained the sum of the A and tt momenta to be equal to the S~ momentum. 
The constraint on the A decay products to give the nominal A mass value 1115. 683±0. 006 
MeV/c^ [5] was also apphed. 

The resulting Att" invariant mass spectrum after the fit is shown in Figure 3 (b). 

4.2 Analysis of the Stt system 

Each reconstructed S~ candidate in the mass range between 1.30 GeV/c^ to 1.34 
GeV/c^ was combined with a pion. 
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The mass spectrum of neutral combinations S~7r"'" is shown in Figure 3 (c); a clear 
5(1530) peak of 820 ± 50 events is observed. The production properties of 5(1530) have 
already been measured by DELPHI in [16]. 

The mass spectrum of combinations 5~7r~ is shown in Figure 3 (d). No significant 
excess is observed. The histogram shows the prediction of the simulation for the 5~7r~ 
spectrum without pentaquarks. To estimate the number of pentaquarks we performed a 
fit of the form (1) to the 5~7r~ mass spectrum, with a Gaussian central value of 1.862 
GeV/c^ and a width of 0.015 GeV/c^ equal to the resolution in this mass region. The 
number of events resulting from the fit is equal to —50 ± 75, dominated by the error from 
the fit itself. The reconstruction efficiency of a possible $(1860) object decaying into 
5~7r~ has been computed from a Monte Carlo generated sample of $(1860) events, to 
be (10.0 ± 0.5)%; the error is dominated by the uncertainties on particle reconstruction 
and identification. This leads to an estimate of the upper limit of the production rate of 
a $(1860) object decaying into 5~7r~ per hadronic Z decay, at 95% CL: 

(A^$(i860)~) X Sr($(1860)— ^ 5-7r-) < 2.9 x 10"^ . 

5 Search for Charmed Pentaquarks in the D*p 

system 

5.1 Event Selection 

After the standard hadronic event selection criteria listed in section 3.1 were applied 
to the data collected in 1994 and 1995, about 2.1 miUion hadronic events remained. 

Events containing the decay chain D*"*" — * D^'X K~7r"'"X were selected as a first step 
of the analysis. The following selection criteria were required to suppress the background: 

• X£'(K7r) > 0.15, where is the energy fraction with respect to the beam energy; 

• in the reconstructed D*' decay, it was required that both the kaon and pion momenta 
were larger than 1 GeV/c, and that the angle between the K and tt momenta were 
smaller than 90° in the D* system; 

• the momentum of the bachelor pion (the soft pion coming from the D* -^Dtt decay) 
had to be between 0.3 GeV/c and 2.5 GeV/c, and the angle between the bachelor tt 
momentum in the rest frame of the reconstructed D° and the momentum of the D° 
candidate had to be smaller than 90°; 

• the decay length of the D° had to be smaller than 2.5 cm, but positive by at least 
three standard deviations; 

• cosOk > —0.9, where cos9k is the angle between the D" flight direction and the K 
direction in the D° rest frame; 

• the invariant mass of the Ktt system had to be between 1.79 GeV/c^ and 1.91 GeV/c^, 
and the mass difference AM = Mk^tt — Mxn was required to be between 
0.1425 GeV/c2 and 0.1485 GeV/c^; 

• the K and n candidates were required to have at least one hit in the VD; 

• the K candidates should not have a positive pion tag. This requirement suppresses 
about 50% of the combinatorial background surviving all other cuts. 
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Figure 4: (a) Invariant K+tt mass, (b) Distribution of AM = Mk-kw — Mx-n-- 

5.2 Analysis of the D*p system 

The Mxn and AM spectra obtained after the cuts hsted above are shown in Figure 4. 
The backgrounds around the very clear D° and D* (corresponding to the decay D* D^n) 
peaks are quite small. 
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Figure 5: Invariant masses (a) M(D*+p) and (b) M(D*+p). 



Figure 5 shows the invariant mass distributions of D*p, for total charge zero (right 
charge for a possible pcntaquark) and total charge 2 (wrong charge) respectively. No 
narrow resonance peak around 3.1 GeV/c^ is seen in Figure 5 (a), which corresponds to 
the right charge. 
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To obtain an upper limit for the production of a possible 0c(31OO)'^ state, a pentaquark 
signal was simulated; the detection efficiency for a 6c(3100)° state decaying into D*~^p 
was estimated to be about 0.8%, taking into account the relevant branching fractions of 
the D* and of the D°. 

The best fit to the mass distribution of right charge pairs for a mass of 3100 MeV/c^ 
and a width corresponding to the experimental resolution, with the same procedure as 
described in the previous sections, gives an excess of 7 ± 4 events. The systematic un- 
certainties, dominated by the uncertainties on particle identification efficiencies, are neg- 
ligible with respect to the error from the fit. The 95% CL upper limit on the average 
production rate, per hadronic Z decay, of a 6c(3100)'^ object decaying into D*~^p, is 

(A^ee(3ioo)o) X Br(ee(3100)° ^ D*+p) < 8.8 x 10"^ . (4) 

6 Conclusions 

A search for pentaquarks in hadronic Z decays was performed, and none of the states 
searched for was found. Upper limits were estabhshed at 95% CL on the average produc- 
tion rates {N) of such particles and their charge-conjugate state per hadronic Z decay: 
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These limits improve previously published results [9]. 

In recent years thermodynamical [17] and phenomenological models [18,19] have ap- 
peared, which successfully describe the overall particle production rates in high energy 
interactions with very few parameters. According to the model by Becattini [17], the 
average production rate for the production of the 0+ at the Z energy should be of 0.007. 
According to the model by Chliapnikov and Uvarov [18], the average production rate is 
expected to be less than 5 xlO~^, if the is dominantly produced from the intermediate 
N*/A* baryon state with the mass of 2.4 GeV/c^ as indicated by the CLAS experiment 
[3]. On the other hand, if the 0"*" production mechanism is similar to the one for ordinary 
baryons produced at LEP, its average production rate should be comparable with that 
of a known resonance, the A(1520), which is observed with an average production rate of 
0.0224 ± 0.0027 per hadronic event [5]. 
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